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Abstract 

We write down the relevant field equations for a stationary axially 
symmetric rigidly rotating dust source in such a way that the general 
solution depends upon the solution of an elliptic equation and upon 
harmonic functions. Starting with the dipole Bonnor solution, we built 
an asymptotically flat solution with two curvature singularities on the 
rotational axis with diverging mass. Apart from the two point singu- 
larities on the axis, the metric is regular everywhere. Finally, we study 
a non-asymptotically flat solution with NUT charge and a massless 
ring singularity, but with a well-defined mass-energy expression. 

PACS numbers: 04.20.-q, 04.20.Jb, 04.40.Nr 

Introduction 



The problem of building a physically admissible metric for an isolated ro- 
tating body is still long an unresolved problem [HE]. In fact, in order to 
obtain a physically reasonable source, many restrictions must be imposed 
(energy conditions, regularity, reasonable equation of state). In particular, 
Einstein's equations for a rotating body with a perfect fluid source seem 
not to be integrable. A remarkable exception is given by dust pressureless 
stationary axially symmetric spacetimes. In a remarkable paper [10], Wini- 
cour showed that Einstein's equations for a stationary axially symmetric 
dust source with differential rotation can be reduced to quadratures. These 
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equations contain as a subclass the van Stockum one [TT] of rigidly rotating 
matter. The first asymptotically flat solution of the van Stockum class can 
be found in [12]. The author in [12] shows that the van Stockum class of 
solutions that are not cylindrically symmetric cannot exsist in the Newto- 
nian theory. The solution in [T2] has a curvature singularity with diverging 
mass. Further, the technique named "displace, cut, reflect" [13J to obtain 
rotating discs immersed in rotating dust must be noted. Unfortunately, this 
method generates distributional exotic matter on the z = plane. 
In this paper, starting from the Lewis [141 [T5] form of the metric, we write 
down the equations for stationary axially symmetric rigidly rotating space- 
times in a co-moving reference frame in such a way that the general solution 
depends upon the solution of an elliptic equation and upon harmonic func- 
tions. The class of solutions contains the van Stockum line element for a 
suitable choice of the harmonic function. 

In this context, starting from the dipole Bonnor solution |12j , we obtain an 
asymptotically flat solution with two curvature singularities on the rotation 
axis and showing simalar properties to the Bonnor solution. 
In section 1 we derive the basic equations. In section 2 we present our 
solution. Section 3 collects some final remarks and conclusions. In the ap- 
pendix we derive a non asymtotically flat solution with a NUT charge and 
a well-defined mass-energy expression. 

1 Basic Equations 

Our starting point is the Lewis [14] line element for a stationary axisym- 
metric space-time: 



ds 2 = e vip ' z) [dp 2 + dz 2 ] + L{p, z)d<j) 2 + 2m(p, z)dtd<j> - f(p, z)dt 2 , (1) 



where x 4 = t is the time coordinate, x 1 = p is the radial coordinate in 
a cylindrical system, x 2 = z is the zenithal coordinate and x 3 = <fi is the 
azimuthal angular coordinate on the plane z = 0. Also, 



Further, the root square of the determinant of the 2-metric spanned by the 
Killing vectors dt,d$ is 



Expression ([3]) characterizes the measure of the area of the orbits of the 
isometry group. In the vacuum, the field equations for ([!]) imply that W(p, z) 



t G (— oo, oo) , p 6 (0, oo) , z G (— oo, oo) , <j> G [0, 2tt). 





(3) 
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is harmonic, i.e. W >a ,a = 0, where subindices denote partial derivative and 
a summation with respect to a = p, z is implicit. Therefore, the function 
W(p, z) can be chosen as a coordinate. Looking for regular solutions on the 
axis, the simplest assumption can be made by setting W = p. In this way, 
the van Stockum line element emerges by taking a dust source. However, 
this is not the most general choice. Thanks to the gauge freedom, we can 
take 

fL + m 2 = p 2 H(p,z), (4) 

where H(p, z) is a sufficiently regular function to be specified by the field 
equations. 

We consider a perfect fluid T^ = (E + P)u^u u + Pg^u, with E being the 
mass-energy density , P the hydrostatic pressure and the 4-velocity of 
the fluid. We consider a co-moving reference frame: 

' ' (5) 
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Denoting with the Ricci tensor, the relevant field equations are 

Rzz Rpp — 0, 

Rpp + Rzz = (P 
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(6) 
(7) 
(8) 

(9) 
(10) 

(11) 



where T = 3P — E. Equation ([9]) involves a second-order partial equation 
for L(p,z), while (fTOl) and (fTTI) give second order equations for m(p,z) and 
f(p,z) respectively. Thanks to fl3J), equations (f9|)- ([TT|) are not independent. 
Therefore, from ([4]), we can express L(p,z) in terms of (f,m,H). Putting 
this expression in ([9]) and using equations (fTU|) and (fTT|) . we obtain the 
following compatibility equation: 



AHPe 1 



H, 



H 2 2 
2H p ' p 



(12) 



In what follows we study dust solutions for which P = 0. Setting p 2 H = F 2 , 
equation ([12]) becomes F >a>a = 0. Thus the compatibility condition for Q- 
(|lip requires that F(p, z) be a harmonic function. Conversely, with F(p, z) 
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no more harmonic, the line element ([I]) is appropriate to describe spacetimes 
with non-vanishing pressure P. For F = p the van Stockum line element is 
regained together with the Papapetrou form of the metric [15] . Equations © 
and ([8]) are linear first-order equations involving v tP and v z and they permit 
us to calculate v >p ,v >z in terms of (f,m,F). By applying the integrability 
condition (v tPjZ = v >ZtP ) for the equations so obtained, we read 

f, z F p = F z f, p . (13) 

Finally, when expressions for v tP ,v tZ are put in ([7]), we obtain 

F,zf,z = -f, P F p . (14) 

Excluding the case F = const (it can be see that this leads to the trivial 
solution E = 0), we have / = const. We naturally choose / = 1. Therefore 
our system of equations is 

F, a , a = 0, (15) 
m~ a - m '" - = 0, (16) 



e v m 2 



E = ^^> ( 17 ) 



[m 2 z F p - m 2 p F tP + 4FF z F ZtP - *FF p F ZtZ - 2m tP m, z F }i 

IFF 2 



,(18) 



[4FF 2 F 2;2 + ±FF p F t2>p - F^ z m 2 z + F )Z m 2 p - 2m, p m, z F p ] 
V ' z ~ IFF 2 kW) 



L + m z = F z . (20) 

First of all, for non-expanding spacetimes, the shear = ^[u^k + u^-i] 
vanishes identically for (fT5]) - ([20|) . and therefore our system of equations de- 
scribes rigidly rotating sources in a co-moving reference frame. When F = p, 
equation (|16p is invariant under the transformation z — » z + a (a a constant) 
and a solution can be expanded as ^m(/),z + aj). Setting F ^= p, if 
F(p, z) , m(p, z) are solutions, then also F(p, z + a),m(p, z + a) are, and thus 
the solutions cannot be expanded. 

Note that we have identified (p, z) with the radial and the zenithal coordinate 
respectively in a cylindrical coordinate system. According to this assump- 
tion, some conditions must be imposed. Firstly, by setting E = 0, the metric 
must reduce to the standard flat expression ds 2 = dp 2 + dz 2 + p 2 d(j) 2 — dt 2 . 
Therefore, \\m.E^F = p. Further, looking for regular spacetimes on the 
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rotation axis , the norm of the space-like Killing vector must be van- 
ishing (except at isolated points) at p = 0, i.e. lim p ^o L = 0. Finally, for 
asymptotically fiat spacetimes, at spatial infinity F(p, z) looks as follows: 
F = p + o{\). 

2 Generating an asymptotically flat solution 

2 

Our starting point is the Bonnor dipole solution 12J F = p,m = — — — 5-. 

(p 2 +z 2 p 

We can obtain a solution of (fT6l) by taking the map p — > F(p, z), z — > G(p, z), 
where F = p (l + j^ij) ,G = z (l - ^Jf^) with c > 0, being b a 
constant. Therefore we get the solution 

F = p( 1 + tJ^), (21) 



m 



(p 2 + z 2 ), 

cp 2 [p 2 + z 2 + bcf 



y / p 2 + z 2 [(p 2 + z 2 ) 2 + 2bcp 2 + b 2 c 2 - 2bcz 2 ] 
c 2 p 2 7(p 2 + z 2 + bcf 



3 ' 
■2 



v = ln(a) + 



E 



8 [(p 2 + z 2 ) 2 + 2/J26C + b 2 c 2 - 2z 2 bc] 
c 2 e- v (3A 



[(p 2 + z 2 ) 2 + 2p 2 6c + b 2 c 2 - 2z 2 bc] 



(3 = a (p 2 + z 2 ) 2 = (p 2 + z 2 ) 2 + c 2 b 2 - 2p 2 cb + 2z 2 cb, 

1 = {p 2 - 8z 2 )(p 2 + z 2 f + 2p A cb + I8p 2 z 2 bc + p 2 c 2 b 2 + 16z 4 6c - 8z 2 c 2 6 2 , 

A = {p 2 + 4z 2 )( /9 2 + z 2 ) 2 - 8z 4 c6 + 4z 2 c 2 6 2 - 6p 2 z 2 cb + p 2 c 2 6 2 + 2p 4 cb. 

Solution ([2T|) is asymptotically flat. Note that the map p — > -F(p, z) , « — > 
(jr(p, z) is not bijective, i.e. is not a diffeomorphism. Concerning the features 
of pij). they depend on the sign of the constant 6. 

For b > 0, apart from p = 0, z = ±Vbc, our solution is regular everywhere. 
At these two points, we have curvature singularities with properties close 
to the p = 0, z = singularity of the dipole Bonnor solution (see [32]). In 
particular, the mass-energy diverges at these points. Otherwise, the energy 
density E(p, z) is integrable. For b = 0, we regain the Bonnor solution. 
Finally, for b < 0, the two point singularities disappear, but emerges a 
curvature ring singularity for z = 0, p = y/\b\c with diverging mass. Inde- 
pendently on the parameter b, at spatial infinity the metric reduces to the 
standard expression in asymptotical cylindrical coordinates, and so also by 
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setting c = (E = 0). Note that, because of the non-invertibility of the map 
between the Bonnor solution and solution (|2ip . the curvature singularity at 
the origin of [12] is shifted in the two curvature singularities of (|2ip (for 
b > 0) on the rotation axis. As a final consideration, it must be noted that 
there exists for (|2ip a finite non singular region about the origin. Thus, 
our solution could be matched, in principle, with some asymptotically flat 
vacuum solution. We do not enter in this discussion, but only mention this 
possibility. 

3 Conclusions and final remarks 

We have studied stationary axially symmetric rigidly rotating dust space- 
times in terms of harmonic functions. In [16] Bonnor found the general 
solution for charged dust with zero Lorentz force in terms of harmonic func- 
tions. However, the use of such kind of functions in [16J is different from 
the one in our paper. In fact, in our paper harmonic functions appear in 
equation (J3|) thanks to the gauge freedom, while in [16] the function F(p, z) 
is chosen to be the cylindrical polar coordinate p. In the Bonnor paper, 
harmonic functions arise in order to obtain the most general solution for 
equation (|16h with F = p. In this case, all the solutions of the equation 
m,ct,a — \ m ,p = are given by taking a generic harmonic function rj(p,z), 
with m = prjp. In a similar way, another harmonic function is introduced 
when charged dust comes in action. Therefore, a direct relation between 
does not exist between the harmonic function F(p, z) and rj(p, z) of [16j . 
Also, the paper [T7] must be noticed in which charged dust solutions are 
given in terms of Bessel functions of first and second kind and hyperbolic 
functions. Also in the paper [T7| the condition F = p is retained. 
In this paper, section two, starting with the dipole Bonnor solution, we 
build a class of asymptotically flat solutions containing the Bonnor one as a 
subclass by a suitable choice of the functions F(p, z), G(p, z). Obviously, it 
must be noted that not all the harmonic functions generate physically sen- 
sible solutions. For a physically sensible solution we mean a regular (apart 
from isolated singularities) asymptotically flat solution. Generally, it is a 
simple matter to verify that, if F(p, z) = p, m(p, z) is a regular differen- 
tiable solution for (fTBj) . then also m(F(p, z),G(p, z)) is a solution for (fTBI) 
with F(p, z) harmonic being G(p, z) the harmonic conjugate to F{p, z) i.e. 
F p = G Z ,F Z = —G p . Further, in order to generate a new asymptotically 
flat and regular solution on the axis starting with a seed solution with these 
two properties, we must build a non-bijective (not a diffeomorphism) map 
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p — ► F(p,z),z — ► G(p,z) such that lim^o = P > limE_>o G = z> an d 
linip^o ^ = linip^o -F 1 = (apart from isolated points) and such that at spa- 
tial infinity the functions F(p, z), G(p, z) look as follows: F = p+o(l) , G = 
z + o(l). Concerning isolated singularities, no general conclusions can be 
made. The functions F(p, z), G(p, z) of section two satisfy all the conditions 
mentioned above. 

Finally, in the appendix we present a non-asymptotically flat solution not 
obtained from a seed solution with the technique discussed above and there- 
fore it represents an ad hoc solution. In particular, it is possible to build ad 
hoc solutions for (fT6|) by setting F(p, z) = p(l + -j^p ) (with c a constant) 
and m(p, z) a homogeneous function such that m iQ;a — ^f- = 0. 



Appendix 

We consider the following solution: 



1 + 



p 2 + z 2 



m 



(22) 



1 + 



16 



ln[p 4 - 2cp z + 2p 2 z 2 + c z + z 4 + 2cz 



--ln[c + p 2 + z 2 } - 2ln[p 2 + z 2 } , E 



c 2 e~ v 



[c + p 2 + z 2 



2 ■ 



with c a constant. When c > 0, the solution (|22p has a curvature ring sin- 
gularity when (p, z) = (-y/c, 0) The axis is regular for z > 0, while it shows a 
conical (no curvature) singularity when z < 0. Further, for z < there is a 
region where closed time-like curves (CTC) appear resulting in a vioaltion of 
causality. However, it is possible to take a simple coordinate transformation 
found in [18], i.e. r = t + 2c0 giving the whole rotational axis free of coni- 
cal singularities. Unfortunately, we are forced to introduce a periodic time 
coordinate r and therefore once again CTC appear. Therefore, the problem 
of violation of causality cannot be avoided with an opportune coordinate 
transformation. 

The spacetime asymptotically reads the expression appropriate for asymp- 
totic NUT metrics [19\ [20] with NUT charge q given by c = 2q. We can 
estimate the mass inside an infinite cylinder of radius R by means of the 
integral 



M(R) 



R 



2n 



dz I dp EFe v c 
lo Jo 



(23) 
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The integral (|23p is well defined everywhere. Thus, for the mass we get the 
formula 



M{R) = 2tt 2 c 2 l + R- VR 2 + 1 • (24) 



Because of the non- asymptotical flatness, the solution (|22[) is not interesting 
in an astrophysical context. However, the natural arena for this solution is in 
the extra relativistic context given by non-Abelian gauge theories or in the 
low energy string theory |211 [22] where, in order to obtain supersymmetries, 
NUT charge comes in action. 
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